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Abstract 

We look for differential equations satisfied by the generalized Jacobi polynomials 
{P"'^''^^'^ i^)} j^^Q which are orthogonal on the interval [—1, 1] with respect to the weight 
function 

where a > -1, P > -1, M >0 and N >0. 

In order to find explicit formulas for the coefficients of these differential equations we 
have to solve systems of equations of the form 

oo 

^ A, {x)D^Pi'^^^^ {x)^Fn{x), 1,2,3,..., 
1=1 

where the coefficients {Ai{x)}°l^ are independent of n. This system of equations has a 
unique solution given by 

This is a consequence of the Jacobi inversion formula 

^ a + P + 2k + l 

f^. (a + /3 + A: + j + l),_j+i ^ 

K—J 

X Pl:r'''-'-'-'^(x)PlC^=^^+'^{x) =5,„J<^,^,J=Q,l,2,..., 
which is proved in this paper. 
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1 Introduction 



In ijT^ T.H. Koornwinder introduced the generalized Jacobi polynomials 
which are orthogonal on the interval [—1,1] with respect to the weight function 

r(« + /3 + 2) , 



oo 





2a+/3+ir(a + l)r(/3 + 1) 



1 - x)"(l + x)f + M5{x + 1) + N6{x - 1), 



where a>— 1,/3>— 1, M>0 and > 0. As a limit case he also found the generalized 
Laguerre polynomials \ L'^'^\x) \ which are orthogonal on the interval [0, oo) with respect 

I J n=0 

to the weight function 

—^x^e-- + M5{x), 
r(a + 1) 

where a > — 1 and M > 0. These generalized Jacobi polynomials and generalized Laguerre 
polynomials are related by the limit 

L^'*^(x)=^limP:'/^AM(^l_|^. 

In Q we proved that for M > the generalized Laguerre polynomials satisfy a unique 
differential equation of the form 

oo 

ai{x)y^^\x) + xy"{x) + (a + 1 — x)y'{x) + ny{x) = 0, 

i=0 

where {ai(x)}^Q are continuous functions on the real line and {aj(x)}^^ are independent of 
the degree n. In [|| H. Bavinck found a new method to obtain the main result of |^]. This 
inversion method was found in a similar way as was done in |Q] in the case of generalizations 
of the Charlier polynomials. See also section 4 for more details. In |^ we used this inversion 
method to find all differential equations of the form 



MY,a^{x)y^'^{x)+NY,h,{x)y^'^{x) + 

2=0 2=0 
OO 

+ MiV^Q(x)y(')(x) + xy"{x) + (a + 1 -x)y'{x) + ny{x) = 0, (1) 

2 = 

where the coefficients {02(2;)}^;^, {bi{x)}'^^ and {cj(x)}^-^ are independent of n and the 
coefficients ao(x), 60 (x) and co(x) are independent of x, satisfied by the Sobolev-type Laguerre 
polynomials \ -L^'*^'^(x) > which are orthogonal with respect to the inner product 

L J n=0 

1 ^00 

<f,9> = TT7—^ I x''e-^f{x)g{x)dx + Mf{0)giO) + Nf{0)g\0), 

00 

n=0 



r(a + 1) Jo 



where a > —1, M > and > 0. These Sobolev-type Laguerre polynomials |L^''^'^(a;)| 
are generalizations of the generalized Laguerre polynomials |-L°'*^(j;)| . In fact we have 

L^'^'O(x) = L^^^'Hx) and L-'O(x) = L(f)(x). 
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In this paper we will prove an inversion formula involving the classical Jacobi polynomials 
which can be used to find differential equations of the form 

OO CX) oo 

Mj2aiix)y^'\x) + Nj2Mx)y^'Hx) + MNY,Ciix)y^'Hx) + 

i=0 i=0 i=0 

+ (1 - x^)y"{x) + [P-a-{a + P + 2)x] y\x) + n(n + a + /3 + l)y{x) = 0, (2) 

where the coefficients {aj(x)}^;^, {hi{x)}'^i and {ci{x)}'^i are independent of n and the 
coefficients ao(x), 6o(x) and cq{x) are independent of x, satisfied by the generalized Jacobi 
polynomials ^. In Q we applied the special case /3 = a of this inversion 



formula to solve the systems of equations obtained in |11]. 

The inversion formula for the Charlier polynomials obtained in [Q] (see also section 4) 
was also used in to find difference operators with Sobolev-type Charlier polynomials as 
eigenfunctions. 

In (I H. Bavinck and H. van Haeringen used similar inversion formulas to find difference 
equations for generalized Meixner polynomials. 

2 The classical Laguerre and Jacobi polynomials 

In this section we list the definitions and some properties of the classical Laguerre and Jacobi 
polynomials which we will use in this paper. For details the reader is referred to |^], |13| and 

Ln {x) \ can be defined by 

J n=0 

n / j_ \ k 

L'nH-) = i:yi)'[l^_iyj^, n = 0,l,2,... (3) 



|17|. 



for all a. Their generating function is given by 

(1 - t)-"-' exp (-^) = £ LtHxr (4) 

n=0 

and for all n € {0, 1,2,...} we have 

D^Li-Hx) = {-iyLi"^\x),z = 0,l,2,...,n, (5) 

where D = — denotes the differentiation operator. The Laguerre polynomials satisfy the 
ax 

linear second order differential equation 

xy"{x) + (a + 1 — x)y'{x) + ny{x) = 0. (6) 

It is well-known that 

^-gM)'(:!:)4°'w. -0.1.2 
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This formula can easily be proved by using definition and changing the order of summation 
as follows 

The classical Jacobi polynomials \ Pn°''^\x)\ can be defined by 

L J 71=0 

P^«(,) = t + ° + (^)\ n = 0.1.2.... (8) 

, (-ir t - ° - '^)'' „ = 0. 1,2, . . . (9) 

= 2-"i:(;;t°)(";^'')(^-i)'(^+ir'. »=o,i.2,... do) 

for all a and /?. For all n G {0, 1, 2, . . .} we have 



D^Pt^Hx)= '''^''^f^^^^ PtV''^'\x),^ = 0,l,2,...,n. (11) 
These Jacobi polynomials satisfy the linear second order differential equation 

(1 - x^)y"{x) + [P-a-{a + f3 + 2)x] y'{x) + n{n + a + P + l)y{x) = 0. (12) 
Further we have for a + /3 + 1 > (compare with [|l5|, page 277, formula (30)) 

2 J 'to {a + P + k + lU^ (X), n-0,1,2,.... (13) 

This formula is much less known than formula (0) for the Laguerre polynomials. However, 
the proof is quite similar. In section 5 we will prove a much more general formula. 
We remark that (|l3|) can be written in a more general form as 

+ 1) ^(a.d), N 

x) 



A (-n)fc(a + P + l)k{a + k + l)n-k{a + /? + 2fc + 1) 



, n = 0,l,2,..., (14) 



to r{a + P + n + k + 2) ^ 

1 /1-x^ 
T{a + P + l) 

which is valid for all a and /3. 

Pn ' (x) > and the Laguerre polynomials < Ln (x) > are 

J n=0 1^ J n=0 

related by the limit 

LH(x)=^lim P("'/^)(1-^). (15) 
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2x 



We remark that if we replace a; by 1 ™ multiply by /3" and let /? tend to infinity 

in the complex plane along the halfline where a + P is real and a + P + 1 > we obtain 
by using (|T5|) and the fact that we have for all n € {0, 1, 2, . . .} 



i-nUa + k + = (-1)'= r ^ ^ j n!, A; = 0, 1, 2, 



,n. 



3 The systems of equations 

Let a > —1. The Sobolev-type Laguerre polynomials |-L°'^'^(x)} can be written as 

L J 71=0 

L^'^'^ix) = AoLI,"\x) + AiDL^^\x) + A2D^4^\x), n = 0, 1,2, . . . , 
where the coefficients and A2 are given by 



An = 1 + M 



n + a\ n{a + 2) — (a + 1) ( n + a 



n — 1 



N 



+ 



(a + !)(« + 3) V^~^/ 

MTV /n + a\/n + a + l^ 



+ 



(a + l)(a + 2)ln-i; V n-2 



\ n I a + 1 \n-ll {a + 1)^ \ n j \ n - 2 



A. 



N (n + a\^ MN fn + a\ln + a + l 



a + l\ n — 1/ (a + 1)^ \ n I \ n — 1 



For details concerning these Sobolev-type Laguerre polynomials and their definition the reader 
is referred to Q and [|l2|. Since the classical Laguerre polynomials ^ satisfy the 

differential equation (^) it is quite reasonable to look for differential equations of the form 
(|l]) for these Sobolev-type Laguerre polynomials ^ in view of this definition 

and the fact that L'^'^'^{x) = L^\x). In it is shown that this leads to eight systems of 
equations for the coefficients {ai(x)}^Q, {^i(2;)}^o ^"^^ {ci(a;)}^o- order to find these 
coefficients we have to solve systems of equations which are of the form 

00 

Y,Ai{x)D'+''L^^\x)=Fn{x), n = k + l,k + 2,k + ?,,..., 

i=l 

where A; S {0, 1, 2, . . .} and the coefficients {Ai{x)}°^^ are independent of n. In it is pointed 
out that this system of equations has a unique solution given by 



A,{x) = {-ir'Y.L. 



(.-x)Fj+k{x), i = 1,2,3, , 



This is an easy consequence of the Laguerre inversion formula 



k=j 



(16) 
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oo 
n=0 



which was found by H. Bavinck in j^]. For more details the reader is referred to and 
See also section 4 of this paper. 

Now we take a > —1 and /3 > —1. The generalized Jacobi polynomials 
can be written as 

pa,/3,M,7V(^) = AoPj^^'^Xx) + [Ai{l - x) - ^2(1 + x)] DP^'^'/^^x), n = 0, 1,2, . . . , 
where the coefficients Aq, Ai and A2 are given by 

fn + p\fn + a + p + /n + a\/n + a + /3 + l\ 



+ MN 



(a + /3 + 2)2 /n + a + /3 + l\^ 



(a + l)(/3 + l)V ^-1 



n + (3\ (n + a + 



^ _ M \ n ) \ n J MN /n + a + /3\/n + a + /9 + l 
^ a + P + l /n + a\ a + l\^n-iy\^ n ^ 



/n + a\ / n + a + /3\ 

iV \^ n y \^ n y MN fn + a + p\fn + a + P + l' 
^2 = 7; — : 7 ^ 1- 



a + P + 1 


ir) 







P + 1\ n — 1 I \ n 



Here we used the same definition as in [14|, but in a slightly different notation. The case 
a + P + 1 = must be understood by continuity. In view of this definition and the fact 
that the classical Jacobi polynomials ^ satisfy the differential equation (^) it 

is quite natural to look for differential equations of the form (|2|) satisfied by these generalized 
Jacobi polynomials as was already pointed out in [l^]. Again this leads to eight systems of 
equations for the coefficients {aj(x)}^Q, {bi{x)}'^Q and {cj(x)},^Q. In order to find these 
coefficients we have to solve systems of equations which are of the form 

00 

J2Mx)D'Pt^\x) = Fnix), n= 1,2,3,..., 

i=l 

where the coefficients {Aiix)}'^^ are independent of n. This system of equations has a unique 
solution given by 

A[x)-2 ^^^^^^^.^^^^_^^P^_j [xjl'.ix), 1-1,2,6,.... 

This is a consequence of the Jacobi inversion formula 

^ g + /3 + 2fc + 1 
^. {a + P + k + j + l\.j+i 
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X P 



{-a-i-l,-l3-i-l) 



i—k 



{a+j,f3+j) 



k-j 



(x) = 6ij, j <i, i,j = 0, 1,2, . . . , (17) 



which win be proved in this paper. Again, the case q + /3 + 1 = must be understood by 

2x 

continuity. We remark that if we replace x by 1 — 3" i'^ (0)' multiply by and let (3 tend 

P 

to infinity along the positive real axis we obtain the Laguerre inversion formula ( |16D by using 

(H). 

In [|l^] we found all differential equations of the form 



M^ai(x)y(*)(x) + (1 - x'^)y"{x) - 2{a + l)xy' (x) + n{n + 2a + l)y{x) = 0, 



(18) 



i=Q 



where {aj(x)}^Q are continuous functions on the real line and {ai{x)}'^^ are independent 
of n, satisfied by the symmetric generalized ultraspherical polynomials 
defined by 



oo 
n=0 



ja,a,M,M 



(x) = CoPt'"\x) - C7ixDP^"'°)(x), n = 0,1,2,, 



where 



Ci 



1 + 



2Mn (n + 2a + l 



a + 1 



n 



n + 2a + l 
n - I 



2M (n + 2a\ 2M'^ ( n + 2a\ ( n + 2a + 



2a + 1 



+ 



n 



a + l\ n — 1 



n 



We remark that these polynomials form a special case {(3 = a and N = M) of the generalized 
Jacobi polynomials |p°'^'^^'^(x)| ^, but the differential equation ([l^) has a very special 
form without a M^-part. This is explained by the fact that 



Co 



1 + 2M 



/n + 2a + 1^ 



\ n — 1 

2M (n + 2a^ 



2a + I 



n 



1 + 2M 



n + 2a + 1 
n - 1 



This implies that the generalized ultraspherical polynomials satisfy the same differential equa- 
tion as the polynomials { Q?!.'^'^^ '^Ux) \ defined by 

L J n=0 



1 + 2M 



+ 2q! + 1 



n — 1 



pt"Hx) + 



2a 



2M (n + 2a\ , 

-( ^ jxI)P^°'-)(x), n = 0,1,2, 



However, this differential equation will appear not to be a special case of the differential 
equation of the form (||) for the generalized Jacobi polynomials, since the MiV-part will not 
vanish if we take (i = a and = M. We aim to give a proof of this in a future publication. 
In 0] we applied the special case /3 = a of the Jacobi inversion formula ([T7[) to solve the 
systems of equations obtained in O]. 
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4 The inversion formulas 

In Q H. Bavinck and R. Koekoek found the following inversion formula involving Cliarlier 
polynomials 

E ctl\-x)C^l^{x) = 6,j, J <i,i, J = 0,1,2,.... (19) 

k=j 

This formula is an easy consequence of the generating function (see for instance [jl3|) 

oo 
n=0 

In fact we have 

oo oo 

1 = e-^'il + t)^e'^*(l + t)-^ = E Ci''\x)t^ 

k=0 m=0 
OD / n \ 

n=0 \k=0 / 

Hence 

f 1, n = 

E^i'^^wC?(-) = 

fc=o [ 0, n = 1,2,3, ... . 

Now ( p!9D easily follows by taking n = i — j and shifting the summation index. This formula 
was also used in to find difference operators with Sobolev-type Charlier polynomials as 
eigenfunctions. In Q a similar formula involving Meixner polynomials was used to find 
difference equations for generalized Meixner polynomials. 

Formula ( [TqI ) can be interpreted as follows. If we define the matrix T = {tij)fj^Q with 
entries 



tij 



0, j > i, 

then this matrix T is a triangular matrix with determinant 1 and the inverse U of this matrix 
is given by = U = {uij)^j^Q with entries 

ct;\-x), j<i 

0, j > i. 

Therefore we call ( p!9| ) an inversion formula. 

In the same way we find by using the generating function (Q) for the Laguerre polynomials 

i 

Y,L^^X{x)l[-_X^\-x) = 5,„ j<i, i,j = 0,1,2,.... (20) 

k=j 

However, this formula cannot be used to solve systems of equations of the form 

oo 

E Aix)D'Ll^^ (x) = Fnix), n = 1,2,3,... 

1=1 
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in view of the parametershift in (|5|) . 

In [HI H. Bavinck used a slightly different method to find the Laguerre inversion formula 
( [16| ) from the generating function (^) for the Laguerre polynomials. In fact we have 

(1 - ty-^-^ = (1 - ty-^-'^ exp ( (1 - exp ^ "^^ 



t-lj ' ' \t-\ 



CO oo 
fc=0 m=0 

n=0 \A:=0 / 

This implies, by comparing the coefficients of t*"-' on both sides, that 



E4"'''^(x)Lt;T ^(-x) = 5.„ J<^, ^,J = 0,1,2,, 



fc=o 

which is equivalent to (jl^). 

Formula (16) can be interpreted as follows. If we define the matrix T = {tij)^j^Q with 
entries 

t'ij = < 

I 0, j> i, 

then this matrix T is a triangular matrix with determinant 1 and the inverse U of this matrix 
is given by = U = {uij)fj^Q with entries 

0, j > i. 

In case of the Jacobi polynomials the above methods seem not to be applicable. In that 
case we have to find the inverse of the matrix T = (tij)"j=Q with entries 

tij = D^P,^'''^\x), i,j = 0,1,2,..., n. 



This matrix T is also triangular and by using (11) the diagonal entries equal 

tu = D^Pr\.)=^^±^^, ^ = 0,1,2,..., n. 

This implies that the determinant of T is nonzero for each n iff — (a + /? + 2) ^ {0, 1, 2, . . .}. 
In that case T is invertible and if the inverse U is given by T"^ = U = {uij)i'j=Q then we 
must have 

1 2* 
Uii = — = 7 ^, i = 0, 1, 2, . . . , n. 

tii {a + f3 + i + iy 

In the next section we will give a proof of the Jacobi inversion formula (0), which is equivalent 
to 

(a + /? + 2j + l)y 
(« + /3 + J + l)m^-^' 

0, j > i. 
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5 Proof of the Jacobi inversion formula 



In this section we will prove that 

^ (« + /3 + 2/c + l)(a + /? + l)fc (-n-a-l-n-P-l) ( . 

,t^o na + P + n + k + 2) ^""'"^^-^ 

' n = 0,1,2,..., (21) 



r(a + /? + l)n! V 2 

which holds for all a and /?. 

Note that (jl|) is a special case of ( pl| ) since 

_„_/3-i)..^ ^ (-n-a)„_fc _, fc (a + fc + l)n_fc 
"-^ ^ ^ (n-A;)! ^ ^ [n-k)\ 



-n)k{a + k + l)n-k, k = 0,1,2, ... ,n 



for ah n G {0,1,2,...}. 

By taking y = a; in (^TJ) we easily obtain 

(a + /3 + 2fc + l)(a + /? + l)fc ^ 
r(a + /3 + n + A; + 2) ^ 

i n = 

r(a + /3 + l)' ^22) 

0, n = 1,2,3,... 

for all a and /3. If we take n = i — j in (p2|) and shift the summation index we find 



^ (g + /3 + 2fc - 2j + l)(a + /? + 1)^,^ ^ 
^. T{a + f3 + i-2j + k + 2) 

For a and /? real with a + /3 + l > —1 we now obtain ( |T7|) by shifting both a and /3 by j. 
Note that ( ]2l| ) for y = — x in a similar way leads to 

^ a + p + 2k + l 



Ua + P + k + j + l)i 



This formula was used in 

In order to prove ( ^l] ) we start with the left-hand side, apply definition (|8|) to pjf^'^\x) 
and definition (^) to i'^.^ a i, n /3 i) change the order of summation to obtain 

^ {a + P + 2k + l)(a + /? + 1)^ p(a,/3)^^^ 

A. r(a + /3 + n + A: + 2) ^^^^n-fc 1?/^ 
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V V Vf ..n-fc (c^ + /3 + 2A: + !)(« + /? + l)fc(a + /3 + fc+l),(a + z + l)fc_, 

r(a + /3 + n + fc + 2) i\ {k - i)\ 

{a + l3 + n + k-j + 2)j{a + k + l)n~k-i ( x-l V ( y-^ V 
j\ {n-k-jy. [ 2 J \ 2 J 

n n n—k 

= EEE(-ir-'x 

i=0 k=i j=0 

{a + P + 2k + l)(a + /? + l)i+k{a + i + l)n-i-j fx-lV fy-l\^ 



T{a + l3 + n + k-j + 2)i\{k-i)\j\{n-k-j)\ \ 2 

n n—i n—i—k 



= (-1) 



n—i—k 



X 



i=0 A;=0 j=0 

(g + /9 + 2i + 2/c + + /3 + l)2i+k{a + ^ + f x - ly ( y - I 



X 



= EE E (-1) 

j=Oj=0 fc=0 



T{a + P + n + i + k - j + 2)i\k\j\{n - i - k - jy. 
X — IV/?/ — IV 



n n—in—i—j 

n—i 



{a + 13 + 2i + 2k + l){a + p + l)2i+k{a + i + l)n-i-j{-n + i + j)k 

X 



= EE(-i) 



r(a + (5 + n + i-j + k + 2)iljlkl{n-i-jy. 
(« + /? + l)2.(a + i + l)n-^-, - 1 V - 1 V ^ 



i=Oi=o ^!i!(n-^-i)! 

^ g r(a + /j + n + i-i + fc + 2)fc! (" + ^ + '^ + '^ + ^^'" = "'^''' 
Now we will show that for all b we have 

In order to prove this we use the well-known Vandermonde summation formula 

-n, b 



2F1 



(c)„ 7^0, n = 0,1,2,..., 



(c) 

which can be written in a more general form as 

{-n)k{b)k ^ (c - b)n „ ^ n . ^ 
t,ric + k)k\ r(c + n)'^ U,l,2,.... 

This formula is valid for all b and c. By using this we find that for all b we have 

^^r{binik + i)k\^^^'^''^ 

_ . V" (-T^)fc(b + l)fc _ , {-n + l)k{b + l)k 



^r(6 + n + A;+l)A;! ^^T{b + n + k + 2) k\ 
no .^^j ^ ^ — = 0, n = 1, 2, 3, . . . , 



r(6 + 2n + 1) T{b + 2n + 1) 
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which proves (|2|). Now we use (p3|) to obtain 



^ {a + P + 2k + l)(a + /? + 1)^ {a,l3), (_„_a-l . s 

r(a + /3 + n + A: + 2) W^n-fc [V) 

{a + P + l)2i fx-iy fy-iy-' a + /3 + 2i + l 



E(-i) 



j=0 



i\ in — i)\ 



T{a + (3 + 2i + 2) 



1 1 " 

r(a + /3 + l)n!^ 



1 



1 f X — y 



r(a + /3 + l)n! V 2 
which proves (|2l|). 



n = 0,1,2, 



6 Some remarks 

Note that we have from definition (^) for the Laguerre polynomials that 

4-"^(^) = (-ir^, n = 0,1,2,.... (24) 
n! 

Hence, the polynomial Ln "'\x) reduces to a monomial of degree n for all n G {0, 1, 2, . . .}. 
Definition ( |lO|) for the Jacobi polynomials leads to 

p<-«(.) = (";'')(i^)".„=o.i.2 

which is also a monomial. However, this monomial might reduce to the zero polynomial. For 
instance, Pn '^\x) equals the zero polynomial for all n G {1, 2, 3, . . .}. 
It is possible to generalize the Laguerre inversion formula ( p!6[ ) to 



f]4"+P")(x)4X'"V^) = ^^^-^^f^, ^ = 0,1,2,..., (25) 
k=o ^■ 

where pn and g„ are arbitrary and even may depend on n. In order to have an inversion 
formula we have to choose pn and g„ such that 

{Pn - gn + 2)n = 0, n = 1, 2, 3, . . . , 

hence 

Pn-Qn^ {-n - 1, -n, . . . , -3, -2}, n = 1, 2, 3, 

Note that the endpoint-cases Pn — Qn = —n — 1 and pn — Qn = —2 correspond to the earlier 
mentioned inversion formulas (|l^) and (^) respectively. 
To prove (B^) we use (H) to obtain 



L(r)(x) = 1 A" [(1 - i)"""' exp (^^) 



, n = 0,l,2,... 

t=o 



11 



where Dt = — denotes differentiation with respect to t. Hence by using Leibniz' rule we find 
at 



k=0 

n 



1 



k=0 



-o-p„-l 



exp 



xt 
t-1 



1 



{n-k)l 



n~k 



(1 - 
{Pn - gn + 2) 



exp 



-xt 



t - 1 



t=0 



-, n = 0,l,2,..., 



which proves (p5|). 

2x 2y 

Further we remark that if we replace x by 1 and y by 1 in (|2TD, multiply by 

f3 _ 
r(a + /3 + l)/?" and let /? tend to infinity in an appropriate way we obtain by using (Qq) 



n—k 



-y) 



{y - xY 



k=0 



, n = 0,l,2,.... 



(26) 



Note that (0) is a special case of ( p6| ) since 

1)(0) = (_i)n-/c (« + fe+^l)n-. ^ (_l)n(_i). f ^ + f ) , = Q, 1, 2, . . . , n 



(n- A;)! 



n — k , 



for all n € {0, 1,2,...}. Moreover, note that ( p6| ) is a special case of the well-known convolu- 
tion formula for the classical Laguerre polynomials 



E Lrix)Lf_,{y) = L(r+^+i)(x + y), n = 0, 1, 2, . . . 
fe=o 

in view of (|2^). By using the technique demonstrated above this convolution formula can be 
proved for all a and (3 which might even depend on n. 
Finally we remark that, by using the fact that 

{h/2)k{h + 2k) = 6(6/2 + l)fc. A; = 0, 1, 2, . . . , 

formula ([2^) can also be obtained by using a summation formula for a terminating well-poised 



hypergeometric series (see for instance formula (III.9) in [16|). 
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